How general are time-local master equations? 
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Time-local master equations are more generally applicable than is often recognised, but at first 
sight it would seem that they can only safely be used in time intervals where the time evolution 
is invertible. Using the Jaynes-Cummings model, we here construct an explicit example where 
two different Hamiltonians, corresponding to two different non-invertible and non-Markovian time 
evolutions, will lead to arbitrarily similar time-local master equations. This illustrates how the time- 
local master equation on its own in this case does not uniquely determine the time evolution. The 
example is nevertheless artificial in the sense that a rapid change in (at least) one of the Hamiltonians 
is needed. The change must also occur at a very specific instance in time. If a Hamiltonian is known 
not to have such very specific behaviour, but is "physically well-behaved" , then one may conjecture 
that a time-local master equation also determines the time evolution when it is not invertible. 

PACS numbers: 03.65 Yz, 42.50 Ct 



I. INTRODUCTION 

In Nature, quantum systems are never truly isolated; 
the dynamics of a quantum system depends, to lesser or 
greater extent, on its interaction with an environment 
For an isolated system, the evolution would be unitary, 
determined by the system Hamiltonian. If the system 
is coupled to an environment, and system and environ- 
ment together undergo unitary evolution, then the re- 
sult is non-unitary evolution of the reduced density ma- 
trix of the system. Tracing over the environment in the 
Schrodingcr equation for system and environment gives 
a first order differential equation for the reduced density 
matrix of the system. In this so-called master equation, 
the environment enters the description only as parame- 
ters, which may be time-dependent, but are nevertheless 
not variables. 

If correlations between the system and the environ- 
ment are short-lived, then the system has no memory of 
its previous evolution, other than what is contained in 
the present state of the system itself. The time evolution 
then has a semi-group structure and the master equation 
can be written in standard Lindblad form 

Hi. Such 

time evolution is said to be Markovian, has been much 
investigated and is well understood. 

Non-Markovian time evolution occurs, somewhat 
loosely speaking, when the future evolution of a quan- 
tum system is affected by its past history. This may 
happen if there is feedback from the environment onto 
the system, which influences the future evolution. How- 
ever, non-Markovian behaviour is also possible when 
the environment is unaffected by the system, and there 
thus can be no feedback in the usual sense (if. Non- 
Markovian processes appear in quantum optics [l|, Q, 
solid state physics [f| and quantum information process- 
ing 0- Even though non-Markovian behaviour appears 
in many branches of physics, there is no definition of 
non-Markovianity that is generally agreed upon. Several 
measures of non-Markovianity have been proposed 0-[Hl 
and compared [lU, [l3| . 



A master equation in Lindblad form is time-local. This 
is a useful way of writing master equations since it, in 
particular, enables the use of quantum trajectory meth- 
ods [l4j|. At first sight, it might seem that time- local 
master equations can only describe Markovian evolution. 
Appearances deceive, however, and they can in fact also 
describe non-Markovian behaviour As long as the 
time evolution is invertible (and in some other special 
cases), it can be uniquely described using a time- local 
master equation Q. If the time evolution is not invert- 
ible, then time-local master equations typically cease to 
uniquely determine the time evolution. One must then 
also be careful when using numerical techniques such as 
quantum trajectory methods. 



In this paper we will give an explicit example where 
two different time evolutions, corresponding to two dif- 
ferent Hamiltonians, both lead to the same time-local 
master equation. The time evolutions are explicitly 
non-Markovian and are constructed using the Jaynes- 
Cummings model [l5j with time-dependent coupling. 
The Hamiltonian for one of the time evolutions changes 
rapidly at the specific time when the time evolution be- 
comes non-invertible. We conjecture that if it can be 
guaranteed that the Hamiltonian does not have such "ar- 
tificial" features, then the time evolution is in princi- 
ple uniquely determined by a time-local master equation 
even in cases when the time evolution is not invertible. 
This would broaden the applicability of time-local master 
equations to encompass an even larger class of physical 
time evolutions. 



In Section UH we review concepts related to master 
equations and time evolution, and the Jaynes-Cumming 
model. We then proceed to construct and investigate 
the two time evolutions with the same time-local master 
equation in Section [IIII We finish with conclusions. 
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II. MASTER EQUATIONS AND TIME 
EVOLUTION 



Master equations describe the evolution of a quantum 
system S coupled to an environment E. Under fairly gen- 
eral conditions, a master equation for the density matrix 
of the system takes the form [l| 

p (t) = -- \H, p (t)J + J Ku,t \p (u)} du, (2.1) 

where % is the system Hamiltonian, and the memory 
kernel fC Ut t is a linear map describing the effects of the 
environment on the system. 

The Born-Markov approximation amounts to taking 
the kernel as IC Ui t [p (u)] ~ ICS (t — u) p (it), and assum- 
ing that the coupling between the system and the en- 
vironment is weak. Then there is no memory, and the 
system dynamics are said to be Markovian. This leads 
to a master equation in Lindblad form 0, Q , 



(2.2) 



Y,'^{LlLkp{t)+p{t)L\L k 

k 



2L k p(t)L{ 



where 7fc > are the decoherence rates, the operators 
L k describe the different decoherence channels, and the 
Hamiltonian term may now include effects arising from 
the environment. 

The master equation (|2.2[) has a time-local form, that 
is, the RHS is a linear operator acting on p(t) and does 
not explicitly involve the past history of the system. Nev- 
ertheless, time-local master equations can describe also 
non-Mar kovian time evolution [l|. A simple motivation 
for this is given in . One uses the fact that any physi- 
cal time evolution in quantum mechanics is described by 
a completely positive map [l?! 

p{t) = Mf>(0)} = E Mt)p(0)A{(t), (2.3) 
k 

where A k (t) are time dependent Kraus operators which 
satisfy ^ k A\{t)A(t) = X, guaranteeing that the map 
is trace-preserving. The solution of a master equation 
should be a valid physical time evolution, and hence it 
can be written as a CP map [24[ ■ We then have 



/C S]t [/5(s)] ds 



(K St t o <fi s ) [p(0)] ds 



(£.,* o S o 4,- 1 ) {p(t)] ds, (2.4) 



where we have omitted the trivially time-local Hamilto- 
nian term. From this we obtain the time-local master 
equation 



where the map At can be identified with J Q K, s .t ° 4> s ° 
(j)^ 1 ds. Therefore, a time- local formulation is in principle 
possible as soon as the time evolution is invertible in the 
time interval considered. This includes non-Markovian 
cases. If the time evolution is not invertible in the time 
interval considered, i.e., if at least two different states 
evolve to the same state at some time t, then the the 
inverse map 4>J X does not exist, and the above argument 
breaks down |25( . 

Non-Markovian systems thus admit a description us- 
ing Lindblad-like master equations of the form (|2.2p , but 
with time-dependent decay rates 7fe(t) and decoherence 
channels L k (t) HQ, which may also depend on the ini- 
tial state of the environment and on the initial state of 
the system itself. In what is sometimes called the time- 
dependent Markovian case, all 7fc(i) > for all times. For 
the truly non-Markovian case, one or more of the 7&(i) 
takes negative values for certain time intervals. Negative 
decay rates may, for example, correspond to processes 
where the system is reversing to its initial state. Sev- 
eral measures have been proposed to characterise non- 
Markovianity [8l4li1| . There arc slight differences between 
these. For example, evolution which according to the 
trace-distance measure @ is "Markovian" may be "non- 
Markovian" in the sense that negative decay rates ap- 
pear in the master equation [l2l |. This is because the 
effect of the "non-Markovian" decay channels, as far as 
trace distances are concerned, may be cancelled out by 
other "Markovian" decay channels with positive decay 
rates (H|. 



III. MASTER EQUATION FOR A TWO-LEVEL 
SYSTEM 

We will now construct a non-invertiblc time evolution 
and the corresponding time-local master equation, with 
the aim of investigating the extent of the validity of this 
master equation. Let us suppose that a two-level quan- 
tum system with excited state |e) and ground state \g) 
evolves according to a CP map p(t) = Ai(t)p(0)A\(t) + 
A2(t)p(Q)A\(t), where the Kraus operators are [lj| 



Ai(t)= \g )(g\ + f(t)\ e)(e\, 
Mt)= v / lHJWl<7>(e| 



(3.1) 



Continuity of p(t) implies that f(t) should be a con- 
tinuous function, and also that \f(t)\ = 1. Since 
Sfe A k (t)Ak(t) = 1 for the map to be trace-preserving 
for any initial state p(0), it must hold that < \ f(t)\ < 1. 
For instance, if p(0) = |e)(e|, then 



p(t) = (l-|/(t)| 2 ) |5><<7| + |/(*)| 2 |e>(e| 
More generally, it holds that 



(3.2) 



p{t) = \f(t)\ 2 p ee (o)\e)(e\ + [p g (0) + ^-\f(t)\ 2 )pee]\g}(g\ 



p(t) = A t [p(t)} , 



(2.5) 



-f(t)p eg (0)\e)(9\+r(t)p ge (0)\g){e\. 



(3.3) 
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If we for simplicity assume that f(t) is real, it is possible 
to show that this time evolution corresponds to a time- 
local Lindblad-like master equation [l6| 



[2a+p(t)a- — a+0-p(i) — p(£)a+aJ\ 



(3.4) 

where <x + = \e)(g\ and <r_ = \g)(e\. (For complex /(f), 
the master equation contains an additional term which 
is unimportant for our purposes.) We can thus identify 
the generally time-dependent decay rate, for transitions 
from \e) to \g), as 7(f) = -2f(t)/f(t). The choice f(t) = 
cxp(— "/t), where 7 > is a constant, corresponds to the 
familiar case of exponential decay with a decay rate 7. 
An I /(f) I which is monotonically decreasing in some other 
way corresponds to a decay rate that varies as a function 
of time. If |/(i)| is increasing, then this corresponds to a 
negative decay rate and non-Markovian behaviour. 

The time evolution becomes non-invertible if f(t) =0 
for certain time(s) [Tq| . Any initial state will evolve 
to \g)(g\ at these times. One such choice is f(t) = 
cos(u>t), which as we shall see corresponds to the Jaynes- 
Cummings model on resonance, describing the interac- 
tion of an atom with a cavity field. We then have j(t) = 
2oj tan(wf), which not only becomes negative, but also di- 
verges when the time evolution becomes non-invertible, 
that is, for wt = (n + 1/2)tt where n 6 Z. One easily 
sees that the solution of the corresponding master equa- 
tion is not unique. Another function f(t) = Acos(u)t), 
where < \A\ < 1 and A is a constant, would also give 
7 (i) = 2wtan(wt) and exactly the same master equation 
(|3.4p . but a different time evolution. Conversely, if we 
know that the initial state of the system is something 
other than |fl l )(5 , |, then it is possible to fix the value of 
\A\. However, if the initial state is |<?)(g|, or if this state 
is encountered during the time evolution, then we can- 
not determine \A\ for the subsequent time evolution. The 
system could start its time evolution with one value of 
A, and continue evolving with another value of A after 
having reached the ground state. Such a time evolution 
would be a solution to the time-local master equation in 
p.4[) . just as valid as a solution with an A that remains 
constant. 

We will proceed to explicitly construct Hamiltonians 
corresponding to two time evolutions that correspond to 
the same time-local master equation for this example. 



energy difference between \g) and |e), and the strength 
of the coupling between the two-level system and the 
cavity field. In the interaction picture, and after mak- 
ing the rotating wave approximation (RWA), the Jaynes- 
Cummings Hamiltonian becomes fl5| 



hAa z — ihVl (a+a — a- a') 



(3.6) 



where 2|A| = luq — uj is the detuning of the system. The 
RWA is a good approximation if 2A <C ojq + w and O <C 
u . 

In order to use dimensionless units, we scale variables 
by the time r as 



(3.7) 

(3.8) 
(3.9) 
(3.10) 
(3.11) 



The time evolution is determined by the time-dependent 
Schrodingcr equation 





t 


t = 


* 




T 


U> Q = 


TOJ , 


UJ = 


TOJ, 


n = 


tSI, 


A = 


rA. 



at 



(3.12) 



In the one-excitation subspace, the basis states for the 
atom and the field are \g, 1) and |e, 0). Within this sub- 
space, the time-dependent state can be expressed as 



|tf(£)> =Ce(t)|e ) 0)+c fl (i)| S ,l), 



(3.13) 



where c e r g \(i) are time dependent coefficients. The 
Schrodingcr equation expressed in matrix form is 



. d fc e (t) 

IT I 

dt\c g (t) 



A -in 

m -a 



Ce{t) 



(3.14) 



For instance, for c g (0) = 1 and c e (0) = the solution is 
O 



!*(*)> = 



sin 



(Sfif) |e,0) 



(3.15) 



A 



(ujfjt) + i^- si 



Rt) 



\9A), 



A. Hamiltonian and time evolution 

The Jaynes-Cummings Hamiltonian for a two-level 
atom coupled to a quantised electromagnetic field is given 
by 

Hjc = hw (h^a + -J +-hLj Q d z —ihQ (a + + a-) (a — a^) , 

(3.5) 

where , a are the creation and annihilation operators 
for the field, uj is the frequency of the field, Hujq is the 



where lor = + n 2 is the scaled Rabi frequency. 

By varying A and n but keeping the Rabi frequency 
ujr = \/A 2 + 51 2 constant, we will now construct two 
different Jaynes-Cummings Hamiltonians that lead to 
different time evolutions but arbitrarily similar master 
equations. Evolution on resonance, with A = 0, gives 
ujr = n. An initial state \e, 0) then gives p ee {t) = 
\c e (t)\ 2 = cos 2 (ujRt),p gg (t) = \c g (t)\ 2 = shi 2 (w fl ,t) and 
Peg = Pge = for the reduced density matrix of the atom. 
This will correspond to the first one of the two time evo- 
lutions (and Hamiltonians) in our example. The second 
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time evolution starts off from the same initial state, also 
on resonance. At a time ti = (n + 1/2)tt/qJr, when the 
state is \g, 1), the Hamiltonian is rapidly switched to off 
resonance with A' ^ 0, also adjusting il' in such a way 
that the Rabi frequency uir remains equal to f2. If the 
switch is instantaneous, the reduced density matrix for 
the atom will after this time evolve with 



Pee{t) 

Pggit) 



A 1 2 



|Ce(*)l 



cos 2 (ujRt), 



I'M 2 



= Pge = 0. 



Q,' 



COS 2 (wRt), 



(3.16) 



If the change is instantaneous and occurs at the right 
time, then the master equations for these two situations 
will be identical. However, if the change is not instanta- 
neous, it will affect the master equation. We now pro- 
ceed to numerically investigate how rapid changes to the 
Hamiltonian affect the time evolution of the atom and 
the corresponding master equation. 



B. Evolution with time-dependent coupling 

In order to change from an on-resonance to an off- 
resonance Hamiltonian with fixed ujr, both the coupling 
il and the detuning A must change. Let f2(i) go from a 
maximum value f2 max to a minimum f2 m in- We choose a 
smooth step function to achieve this, 



{1 -tanh [k (t-tj)]} + fi r 



(3.17) 

where k = k/r controls the rate of change and ti is 
the point at which tt(t) takes its average value (ri m ax + 
^min)/2, see Fig. [TJ The detuning changes as A(t) = 
\J uj r — ri 2 (t) , with ujr kept constant. Changing the cou- 
pling strength may be feasible in cavity QED [TU, by 
changing the position of an atom in a laser cavity, or by 
changing the laser field strength. In circuit QED, it is 
possible to realize tunable resonators [ljl [2(J . 

The numerical solution of (|3.14[) with fl(t) given in 



(|3.17p . where k 

^ma.x — 0.3, f! m j 



= 1.6, ti chosen so that c g (t.j) = 0, 
= 0.2, and c g (0) = and c e (0) = 1, 



is plotted in Fig. O It is seen that this produces the 
desired time evolution. If the change in the coupling is 
not instantaneous, then this can be seen in the resulting 
time evolution. Fig.[3]shows p ee {t) for different values for 
k. As the change becomes more rapid, the resulting dy- 
namics approach those for an instantaneous change, see 
Fig. [3] As k decreases, and the change becomes faster, 
the amplitude of the oscillations after the change in the 
Hamiltonian decreases. 

Fourier analysis of c e (t) and c g (t), see Fig.|U confirms 
that the frequency of the oscillations in p ee stays con- 
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FIG. 1: (Color online) Scaled coupling strengths Q = Qt for 
different rates of change, varying according to Eq. (|3.17|) . with 
U = U/t — 0. Here, r2 max = Tf2 max is 0.3 and 0, m i n = rSl m i n 
is 0.2. The scaled values of k — rk are 0.1 for the dotted line 
(blue), 0.5 for the dashed line (red) and 1 for the solid line 
(gold). 
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FIG. 2: (Color online) Time evolution of the reduced den- 
sity matrix elements of the atom for changing detuning A(t) 
and coupling Q(t), with constant ujr. The dotted line (blue) 
shows p gg (t) and the dashed line (red) p ee (t). The parameter 
controlling the rate of change of the coupling Q(t) is given 
by k — 1.6. The solid line (gold) shows how the coupling is 
changing. 



stant to a good approximation. We have also numeri- 
cally confirmed that the decay rate in the corresponding 
time-local master equation is practically indistinguish- 
able from the decay rate for the case when the Hamil- 
tonian is not changing in time. We therefore have an 
explicit example where two different time evolutions cor- 
respond to the same time-local master equation. The 
procedure is nevertheless somewhat artificial in that it 
involves a rapid, ideally instantaneous, change in one of 
the Hamiltonians at a very specific time in order to get 
the desired time evolution and an unchanged decay rate. 
If the coupling and detuning do not change at the right 
time or rapidly enough, then p ee does not reach zero ex- 
actly, and then the time evolution will be invertible. Even 
if p ee does reach zero, if the change in the Hamiltonian is 
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FIG. 3: (Color online) Time evolution of p se (t) for different 
values of k. The dotted curve (gold) corresponds to k — 0.5, 
in which case the system never reaches the ground state. As 
will be seen later, this affects the decay rate. The dashed 
curve (red) corresponds to k = 1.0 and the dot-dashed (blue) 
one to k = 1.6. The solid line (green) corresponds to k — > oo, 
that is, an instant change. 
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FIG. 4: (Color online). Fourier analysis of c e (i), with k = 1.6, 
f^max = 0.3 and Q m i n = 0.2. When k —¥ oo, c e (t) and c g (t) 
are piecewise defined functions, with c e (t) — cos(ijfit) if t < U 
and c e (t) = Acos((jjRt) if t > ti, with < A < 1, and c g (t) = 
\J\ — c e (t)| 2 (up to a phase sign). Fourier analysis of c e (t) 
when k — y oo gives a delta-function peak at lor = tujr — 0.3, 
together with a term proportional to 1/(t[cj — Ur\). 

not instantaneous but occurs over a short but finite time 
period, then the master equation will not remain identi- 
cal to the case where the Hamiltonian does not change, 
remaining on resonance. 



C. How rapidly should the coupling change? 

We now discuss a practical estimate of how quickly 
the coupling should change for the time evolution to be 
experimentally indistinguishable from the case in which a 
change is made simultaneously. If a Hamiltonian changes 




FIG. 5: (Color online) This figure shows the numerically ob- 
tained decay rate, the time evolution of the reduced density 
matrix elements and the change in the coupling for our exam- 
ple with the Jaynes-Cummings model. The dotted line (blue) 
and dashed line (red) correspond to p gg (t) and p ee (t) respec- 
tively. The dot-dashed line (green) represents the tangent-like 
decay rate and the solid line (gold) the time-dependent cou- 
pling with Sl max = 0.3, Qmin = 0.2 and k = 1.6. 

from H(to) = %o to H(ti) = Hi between the times t 
and t\, then the probability that the state of the system 
remains unchanged from its initial state is given by 
1 — £, with 

£ = ^ ((Vo|W 2 |Vo) - (V>o|W|'0o) 2 ) , (3.18) 
where T = t\ — to and 

H=^[ 1 H(t)dt. (3.19) 

1 Jto 

If £ <C 1, then the system evolves diabatically, that is, 
the change in the Hamiltonian is so rapid that the sys- 
tem does not have time to adjust. It is straightforward 
to show that this leads to k 3> (f2 ma x + ^min)/2- Ex- 
perimentally, it may be challenging to implement a rapid 
and precise enough change in the coupling strength and 
detuning. 

D. Decay rate in the time-local master equation 

From the numerical solution of the Schrodinger equa- 
tion for the case where the coupling changes, we can also 
compute the decay rate in the corresponding time-local 
master equation. The result is shown in Figs. [5] and |5] As 
expected, the decay rate as a function of time approaches 
the functional form expected for an instantaneous change 
7(<) = 2LUfj teaiujftt, as the change becomes more rapid. 
This means that by looking only at the decay rate, we 
cannot distinguish between a Hamiltonian that does not 
change in time, and a Hamiltonian with a rapid enough 
change at the right time. From Fig. [SJ we also see that 
the decay rate indeed diverges when the time evolution 
is non-invertible, i.e. when the system is in the ground 
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FIG. 6: (Color online) Decay rates for different values of k. 
The dotted curve (gold) corresponds to k = 0.5, the dashed 
curve (red) corresponds to k = 1.0 and the dot-dashed line 
(blue) to k — 1.6. The solid line (green) corresponds to k — > 
oo. 

level. Also, a negative decay rate corresponds to "re- 
coherence" of the system and non-Markovian evolution. 

IV. CONCLUSIONS 

We have investigated time-local master equations 
through an example involving a two-level system. We 
constructed two different system-environment Hamilto- 
nians, corresponding to two different time evolutions for 
the system, which nevertheless both give the same time- 
local master equation. This explicitly shows that the 
time-local master equation on its own is not enough to 
solve for the time evolution, if the time evolution is not 
invertiblc. 

The example is nevertheless somewhat artificial since it 
involves rapid changes in one of the system-environment 
Hamiltonians at a specific instant in time. If any Hamil- 
tonian is guaranteed to be "physically well-behaved", 
for example, that it is continuous and does not change 
too fast (or more precisely, that its matrix elements are 



Lipshitz-continuous) , one may conjecture that the cor- 
responding time-local master equation docs determine 
the time evolution, at least in principle, even when the 
time evolution is not invertible. This would in other 
words mean that two "physically well-behaved" system- 
environment Hamiltonians, corresponding to different 
time evolutions for a quantum system when its environ- 
ment is traced out, cannot both lead to the same master 
equation for the system. Equivalently, this would mean 
that if a master equation does not have a unique solu- 
tion, then at least one of the solutions corresponds to an 
"ill-behaved" system-environment Hamiltonian involving 
rapid changes. 

Now, as our example shows, even a "well-behaved" 
Hamiltonian, such as the one in the Jaynes-Cummings 
model on resonance, may result in divergencies in the 
decay rate in the corresponding master equation. The 
usual theorems related to the existence and uniqueness 
of solutions of differential equations are of little help in 
proving our conjecture. The Picard-Lindclof theorem, 
for example, just tells us that the solution of such a mas- 
ter equation, with diverging decay rates, is not unique. 
If the decay rates in a Lindblad-like master equation do 
not diverge, then its solution would be unique - but this 
simply corresponds to the case where the time evolution 
is always invertible. When the time evolution is not in- 
vertible, then decay rates will inevitably diverge at these 
times, and as already stated, this can and does happen 
even for very well-behaved system-environment Hamilto- 
nians. 

If the time evolution is not invertible, then even if our 
conjecture holds true, and if the Hamiltonian is well- 
behaved enough for the time evolution to be uniquely 
defined by the time-local in a formal sense, one would 
still need to take care when numerically solving a time- 
local master equation. This is because the diverging de- 
cay rates may lead to instabilities in numerical calcula- 
tions. Nonetheless, our results would generally support 
the view that time local master equations are applicable 
to a wider class of problems than one might expect on 
first inspection. 
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